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ABSTRACT 


Properties  of  two  alternative  procedures  to  the  Jack¬ 
knife  Point  and  Confidence  Interval  Estimation  Procedure  of 
Gaver  and  Chu  have  been  studied.  They  are  called  the  Log- 
Noraal  Likelihood  Procedure  (LNLJ)  and  the  Moment  Procedure 
(HP).  These  two  procedures  were  investigated  and  compared 
with  the  Jackknife  Point  and  Confidence  Interval  Availabil¬ 
ity  Estination  Procedure.  Numerical  results  from  simula¬ 
tions  are  presented  in  this  report. 
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I.  I NTRO  DOCTIOH 


A.  OVERVIEW 

Availability  is  the  measure  of  equipment  effectiveness 
that  relates  reliability  and  aaintainability  to  operational 
readiness.  In  some  cases,  availability  and  operational 
readiness  have  been  considered  the  sane.  In  general,  these 
are  all  requirements  which  nust  be  considered  during  the 
design  and  operation  of  a  system,  and,  thus,  must  be  quanti¬ 
tatively  evaluated. 

There  are  a  variety  of  ways  of  expressing  availability. 
In  general,  availability  relates  up  time  (reliability- 
related)  to  down  time  (maintainability-  related)  ,  and  it  may 
be  defined  as  the  ratio  between  the  time  the  system  is  capa¬ 
ble  of  performing  its  mission  to  the  total  time  the  system 
is  in  operational  demand.  Alternatively,  it  is  a  aeasure  of 
the  probability  that  a  single  system  is  "up"  or  available  at 
a  possibly  '’random"  point  in  time  when  its  services  are 
needed. 

The  three  expressions  of  availability  of  greatest  con- 


See 


Rise  and  B  jorklund[  Ref.  1].  In  this  study,  we  will  be  con¬ 
cerned  with  only  the  inherent  availability.  Often,  inherent 
availability,  which  is  a  hardware  oriented  aeasure,  is  spe¬ 
cified  and  required  within  a  maintainability  contract 
requirement. 

A  matter  of  primary  concern  is  that  of  estimating  avail¬ 
ability  from  observations  on  system  up  times  and  down  times. 
The  estimating  process  should  furnish  both  a  point,  or  sin¬ 
gle  number,  estimate  and  also  some  aeasure  of  the  stability 
of  the  estimate  such  as  a  standard  error,  or  confidence  lim¬ 
its.  The  latter  problem  is  more  difficult  than  the  former. 

In  the  paper  by  Gaver  and  Chu[Ref.  2],  it  has  been 
demonstrated  that  the  jackknife  technique  can  be  useful  for 
confidence  interval  estimation  of  system  availability .For 
instance,  let 

2(0) 

A  - -  (1.1) 

2  (0)  *2(D) 

be  the  long  run  system  availability  of  a  single  unit  system 
that  changes  from  up  to  down  states  in  accordance  with  a 
two-state  semi-Harkov  (or  more  general)  process  for  which  up 
[down]  state  expected  duration  is  2(0)[2(D)].  Then  it  is 


9 


possible  to  successfully  jackknife  the  naive  point  estimator 
that  involves  replacing  expectations  by  the  corresponding 
sample  means, 

0 

-  (1.2) 

o*3 

after  initial  logistic  transformation.  There  are,  however, 
particular  parametric  families  of  distributions,  popular  as 
models  for  summarizing  up  and  down  time  data,  for  which  the 
most  statistically  efficient  estimators  of  E  (0)  and  E  (D)  are 
not  the  simple  sample  arithmetic  means. In  particular,  this 
is  true  for, 

(1) .The  log-normal  distribution  with  unknown  (to  be 

estimated)  parameters;  often  this  is  used  as  a  model 
for  down  or  repair  times.  In  particular,  the  log 
normal  distribution  has  enjoyed  considerable  popu¬ 
larity  for  representing  electronic  system  repair 
times. See  Kline  and  Almog[Ref.  3], 

(2) . the  general  gamma  distribution  with  unknown  shape 

and  scale  parameters,  and 

(3) . the  general  Weibull  distribution  with  unknown  shape 

and  scale  parameters. 
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B.  PURPOSE  AND  APPROACH 


1 .  Objectives 

Availability  estimation  by  use  of  the  jackknife  has 
already  been  studied  by  Saver  and  Chu.This  thesis  presents 
two  procedures  alternative  to  the  above.  These  may  be 
called  the  Log-normal  Likelihood  jackknife  (LNLJ)  procedure 
and  the  Moments  procedure  (MP) .  In  brief,  the  jackknife 
method  has  the  capacity  to  reduce  the  statistical  bias  of 
estimates  of  such  quantities  as  system  availability,  and 
also,  and  more  importantly,  to  furnish  confidence  limits 
that  behave  in  a  satisfactory  manner  (economically  enclose 
the  true  availability)  despite  the  fact  that  underlying  dis¬ 
tributions  are  unknown.  See  Miller[Ref.  4]  for  a  review  of 
much  of  the  literature  on  jackknifing. 

The  present  estimation  procedure  makes  use  of  the 
specific  assumption  of  a  log  normal  model  for  repair  times 
to  compute  the  maximum  likelihood  estimate  (ale)  of  avail¬ 
ability.  Properties  of  the  procedure  have  been  studied  by 
Monte  Carlo  simulation.  A  number  of  such  simulation  results 
are  presented  in  this  thesis.  Also  presented  are  compari¬ 
sons  with  the  direct  method  of  Gaver  and  Chu,  and  with  a 
computationally  simple  method  based  on  moment  estimation. 


■ 

s 

Attention  also  has  been  paid  to  study  of  the  sensitivity  of 
the  LNLJ  procedure  to  specification  error.  This  refers  to 
the  error  introduced  by  assumption  of  a  log-normal  distribu- 
[  tional  model  and  the  corresponding  maximum  likelihood  esti- 

i 

[  mates  when,  in  fact,  another  distribution  governs  the  obser- 

\  vations.  It  is  found  that  LSLJ  method,  which  is 

» 

|  theoretically  most  efficient  in  large  samples  if  the  model 

t 

•  (likelihood)  is  correct,  can  sometimes  Droduce  consistently 

fr 

biased  results  when  the  log  normal  model  does  not  apply. 
Since  it  is  nearly  impossible  to  assure  log  normality  from 

'  * 

small  samples  of  data,  this  finding  suggests  that  caution  is 
[  in  order. 


2.  122  normal  Likelihood  Sioc^dure  Poi  2  Sail 


Furthermore,  the  maximum  likelihood  estimates 
(mle)  of  >,id  and  are 

=  x  sl/nSln^)  (1.2.4) 

a2  =1/n  Z  (x.-x)  *  ,  (1.2.5) 

a  i 

so  the  ale  of  E[D]  is,  by  invariance, 

EC  D  ]  *  exp  (wd  ♦  1/2*  7|)  (1.2.6) 

it  is  essentially  the  latter  expression  that 
will  be  used  in  the  estimate  for  availability  to 
replace  the  simple  first  moment  estimate  E(D)=d. 
In  fact,  <?d  has  been  replaced  by  its  unbiased 
version  (n  replaced  by  n-1). 

(c)  .If  the  above  holds,  it  is  advantageous  to  trans¬ 
form  first  (See  Hosteller  and  Tukey[3ef.  5])  rhe 
estimated  availability 
2(0) 

- - ; 

E  (0)  *E  (D) 

the  log-logistic  transform  of  availability  is 
In  (A  /  1  -  A)  =  lnT(u)  -  (^  *1/2*  a  d*)  (1.2.7) 

Jackknifing  will  be  carried  out  using  the  sta¬ 
tistic 

L  *  In  (A  /  1-A)  =  ln(u)-(x*1/2*s  *)  (1.2.9) 

uN  A 
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(i) , Recompute  L  reoeatedly,  leaving  out 
LN 

sively  the  sample  pairs  lu  ,&  ) ,  (u^. 


<W 


j-1 


ii  =  1/n-1.(I  \*Zh  ) 

J  i  =  l 

x.jVn-Lf  g  x..  J  x.  } 

'•V;,  n 

sx-Tt/n”2#?X(x i“  x.)2+.I(xi“  XJ 


succes- 


(1.2.9) 

(1.2.10) 
(1.1.11) 


L,N  ln(tt)-(xj+  1/2*sx^  (1.2.12) 

where  j*1,2, . ,  n 

(e)  .Compute  the  pseudovalues  as  follows: 

l"1'2 . C  <'-2-'3> 

Recall  that  L,  =L  is  the  result  of  the  coa- 
LN  LN,alf 

puting  the  quantity  to  be  jackknifed,  leaving 
out  none  of  data. 

(f)  .Compute  the  mean  and  variance  of  the  pseudova¬ 

lues. 


P  *  Vn»  f  PTM  . 

LN  *-  LN  ,  3 

J-‘  ft  _ 

S  2=  1/n- 1*  T  (P  -  P  ) 2 

LN  j~  LN , j  LN 


(1.  2.  14) 
(1.2.  15) 


(g)  .The  jackknifed  point  estimate  of  the  availabil¬ 

ity  is  now, 

A  =  exp  (P^)  /  1+ex?  (P^  )  (1.2.16) 

(h)  . Symmetric  two  sided  confidence  limits  at  confi¬ 

dence  level  (1-a)^  are  derived  as  follows; 

CL  *  P  -  t  (r. -  1 )  SQ ST  (1/n. S  2)  (1.2.17) 

LN  1-V2  LN 
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(1.  2. 18) 


CaLN=  tl-a/2(n-1)  SQBT(1/n.S J) 

where  2  )  -s  the  P”  * /2)  100%  quantile  of 

Student's-t  with  n- 1  degrees  of  freedom.  Then, 
expCCL^)  exptca^,) 


-  <  A  <  — - -  (1.2.19) 

1  ♦  exp  (CL^  )  1  ♦  exp  (CO^) 

with  confidence  approximately  (1-a  )  100%.  Note 

that  the  confidence  limits  are  nearly  symmetric 

around  In  (E[  0  ]/E[D  ])  ,  and  not  around  k. 

(i) .One  sided  confidence  limits  at  confidence  level 

(1- a  )100%  are  derived 

CLlns  tl-0  <n_1)  SaRT(1/n*SJ)  (1.2.20) 

COms  Pln+  *l-a  SQRKI/n-SJ)  (1.2.21) 

So,  one-sided  upper  confidence  limit  is 
exp  (CO  ^) 


A  < 


(1.2.22) 


1  *  exP<C0^) 

and  lower  confidence  limit  is 
exptCL^) 


A  <  -  (1.2.23) 

1  ♦  expfCt^) 

Note  that  if  up  times  are  exponentially  distributed 
then  u  is  actually  the  ale  of  2[0],  so  under  the  model 
assumptions  the  (transformed)  ale  of  availability  is  being 


5 


jackknifed;  this  procedure  has  been  validated  by  Reeds 
[Ref.  6]  for  large  samples.  Simulation  studies  are 
nevertheless  essential  for  investigating  properties  of  the 
jackknife  or  alternative  techniques  for  small  samples,  and 
for  investigating  sensitivity  to  model  assumption. 

3.  Moment  Procedure  Fjr  |  Siaaie  Unit 

(a) .  Data  given  as  follows 

«1»u2» . *un  (up  times) 

d^  , d^  . ,  d  (down  times) 

Y1  =lnuL  ,Y2  =lnu2  , - ,yn*lnun 

*1  =lnd1  ,x2  *lnd2  #-.*  ..x^lnt^ 

(b) . Assume  up  and  down  times  are  log-normal. Then, 


ECO]*  exp(vu+  i/2*ou2) 

(1.3.1) 

2[D]=  exp(ud*  1/2*ad2) 

(1.  3.2) 

where. 

uu=  E[ InU  ] 

(1.3.3) 

u  *  E[ InD ] 

a 

(1.  3.4) 

a 2  =  7ar[ln0] 

U 

(1.3.5) 

7arClnD] 

(1.  3. 6) 

Furthermore,  the  mle  of  u  »  :*»  j  , 

u  u  d 

and  j 2  are 
d 

\*  Y=  Vn«  £  Yi 

«  A 

(1.3.7) 

au3  Vr.*  £  (yi-y) 2 
»•!  . 

(1. 3.8) 

V  *  -Vn.  £  xL 

(1.3.9) 

i»« 

16 


(c> 


n 


a2  =  1/n*  £  (X.  -X)2  (1.3.10) 

so  the  mle  of  E[0]  and  E[D]  are,  by  invariance, 
5[tJ]=  exp(tu*  1/2* \2)  (1.3.11) 

2[D]=  exp(~d*  1/2*  1d2)  (1.3.12) 

If  above  holds,  the  logistic  transform  of  avail¬ 
ability  is. 


In  (A  /  ( 1-A)  )  =  (^1/2  au*)-(  yd*1/2  od*)  (1.3.  13) 
estimation  will  be  carried  out  using  the  statis¬ 
tic 

P»  In  (A  /  (1-A)  )  =  (y*1/2*Sy2)-(x  +  1/2»sx2)  ( 1.  3. 1U) 

(d)  Simply  the  point  estimate  of  availability  is 

now, 

A=  exp(P)  /  1+exp (P)  (1.3.15) 

(e)  Variance  of  the  statistic  ?  can  be  calculated  as 
follows, 

Var  (P)  =  (SS)  2=  1/n»sy*  ♦  1/4*Var[sy2)  ♦  1/n*sx*  ♦ 
1/««VarCsx2]  (1.3.16) 

where, 

7arCsy2]*  1/n  [  ti4y- (sy2)  2  ]  (1.3.17) 

Var[sx2]=  l/n  CM4x-(sx2)2]  (1.3.13) 

(f)  Symmetric  two  sided  confidence  limits  at  confi¬ 
dence  level  (1-a)%  are  derived  as  follows. 

Lx  p  -  zW2(SE)2  (1.3.19) 
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(1.3.20) 


0*  P  ♦  *L-a/2<SB)» 
where  Zi-a/2-s  the  t1”  »/2)100X  quantile  of  the 
standard  noraal  distribution.  Then, 

(exp  (L)  /  Hexp(L)}  <  &  <  (axp(0)  /  1*exp(0)} 

(1.3.21) 
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II .  STS2i A2IQN  PROCEDURE 

A  simulation  procedure  has  been  used  to  coapare  the 
operating  characteristics  nf  LNLJ  methodology  introduced 
here  with  other  approaches.  Specifically,  siaulation  has 
been  used  to  compute 

(a)  the  actual  coverage  of  the  true  availability  figure.  A, 
by  the  confidence  intervals  given  by  the  procedures 
under  study,  when  tha>  uv;e*.t4.l  coverage  is  a  )  1 00%, 

(b)  measures  of  confident  v’"srval  expected  width  and  vari¬ 
ance  width, 

(c)  estimate  of  the  expected  point  availability  estimated  by 
the  procedure  under  study. 

The  siaulation  program  was  written  in  FORTRAN  IV,  and 
the  siaulation  have  been  carried  out  on  the  I3M  3033  at  the 
Naval  Postgraduate  School.  A  aore  detailed  description  of 
the  prograa  and  its  major  subroutines  is  given  in  Appendix 
A. 

An  outline  of  the  siaulation  procedure  now  follows; 

(A)  A  simulated  sample  of  up-time  durations,  (ui#  i  =  l, 
2,...,  n)  and  down-time  durations,  (di ,  i=1,  2,...,n)  were 


1 
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times.  The  Naval  Postgraduate  School  LLH ANDOM  package  was 
used,  along  with  the  International  Mathematical  and  Statis¬ 
tical  Library  (IMSL)  random  number  generator.  For  all  cases 
considered  (exceptions  noted),  expected  up  and  down  times 
were  E[0]=73.036  and  E[  D  ]=3 . 844  23 1  for  the  population  sam¬ 
pled,  so  long-run  availability  A=0.95.  Sample  sizes  of  n=15 
and  n=25  were  simulated,  and  provided  the  bases  for  point 
estimators  and  for  confidence  intervals. 

(B)  Here  are  the  different  distributional  situations 
that  were  sampled.  A  total  of  1000  replications  was  used  to 
evaluate  each  procedure  in  each  distributional  situation. 

(1)  Op-times  come  independently  from  the  iid  Expon(x) 
distributions.  Down-times  come  independently  from 
the  log-normal  (y, a2)  distribution. 

(2)  Op-times  are  iid  Expon(x).  Down  times  are  iid 
Expon  (u)  . 

(3)  Op-times  are  iid  Expon(X).  Down-times  are  iid 
Gamma  (e,k),  being  a  scale  and  k  a  shape  parameter, 
so  arranged  that  E[  D  ]=  and  7az[D]-(fJc'^  )-*.  A 
Gamma  with  k>1,  e.g.  k=2,  represents  iata  tha'*:  is 
more  tightly  grouped  around  ins  mean  than  is  true  of 
exponentially  distributed  data;  it  roughly  resembles 
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data  that  is  log-normal.  A  gamma  with  k<l,  e.g. 
k=1/2,  represents  data  of  relatively  extreme  disper¬ 
sion  and  positive  skewness  as  compared  to  exponen¬ 
tial  data  having  a  very  long  right  tail,  compensated 
for  by  a  high  density  near  zero.  A  gamma  distribu¬ 
tion  with  integer  k  (k= 1,2, 3, . . . . )  is  often  called 
Erlang;  realizations  are  easily  simulated  by  summing 
k  independent  Exponential  realizations. 

(4)  Op-times  are  iid  Expon  (A).  Down-times  are  iid 
long- tailed  h  distribution; 

d=E[  D  ]*w*exp(h*w)  •  (1-h)  2  ,  0<h<1, 

where  w  comes  from  Expon(l).  This  long-tailed  h 
distribution  possesses  no  closed-form  representa¬ 
tion.  However,  it  can  be  shown  to  have  characteris¬ 
tics  similar  to  those  of  the  Exponential  for  small 
values  of  h,  but  to  have  a  systematically  longer 
right  tail  than  the  exponential. 

(5)  Op  times  are  iid  long-tailed  h  distribution.  Down 
times  are  iid  Expon  (u)  . 

(6)  Op  times  are  iid  Expon  (A).  Down  times  are  iid  Wei- 
bull  (9, k),  being  a  scale  parameter  and  k  a  shape 


parameter,  so  arranged  that  S[ D  ]=?:/*•  r(i/k  ♦  1)  and 


Var£ D ]*92/k«[  r  (2/k  ♦  1)  -  r*(1A  ♦  1)  ].  A  tfeibull 
with  k>1,  e.g.  k=2,  represents  data  that  is  more 
tightly  grouped  around  its  aean  than  is  true  of 
exponentially  distributed  data;  it  quantitatively 
resenbles  data  that  is  log-normal. A  Weibull  with 
k<1,  e.g.  k*1/2,  represents  data  of  relatively 
extreme  dispersion  and  positive  skewness  as  compared 
to  exponential  data  a  very  long  right  tail,  compen¬ 
sated  for  by  a  high  density  near  zero. 

(7)  Op  times  are  iid  Expon(^)  .  Down  times  are  iid 

long-tailed  log-normal  h  distribution.  See  Appendix 
B. 
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III.  ANALYSIS 


The  methods  described  were  simulated  for  various  distri¬ 
butional  assumptions  which  were  defined  in  the  Chapter  2. 
Simulation  results  for  each  method  and  set  of  distributional 
assumptions  are  shown  in  Table  1  and  Table  2. 

In  general,  HP  has  very  high  coverage  factor  and  satis¬ 
factory  point  availability  estimation.  However,  both  average 
and  variance  of  confidence  length  are  higher  than  JK  and 
LNLJ  procedures. In  addition,  no  consistent  superiority  of 
LNLJ  procedure  over  JK  procedure  has  been  noted. 

Especially  under  the  Exponential  up  and  Gamma  (k=l/2) 
down,  and  Exponential  up  and  Neibull  (k=1/2)  down  times 
assumptions,  a  consistently  low  coverage  factor  and  point 
availability  estimation,  and  a  consistently  high  average 
confidence  length  and  variation  have  been  obtained  from  LNLJ 
procedure.  In  these  cases,  distributional  characteristics 
of  down  times  after  log-transform  provide  highly  left  skewed 
shape  and  large  variance,  and,  also  these  features  reflect 
in  the  pseudovalues  which,  in  turn,  causes  left  skewed  shape 
rather  than  a  symmetric  shape  for  the  pseudovalue 
distribution. 


L.'«  •:  — Hviiinrfc'Mfc 
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In  order  to  compensate  for  this  situation,  the  3iveight 


procedure,  see[Ref.  7],  and  the  Winsorizing  procedure,  see 
[Ref.  ‘ ],  were  implemented  directly  on  the  pseudovalues. 
Results  are  shown  in  Tables  3,  4,  5,  6,  7,  8.  These  proce¬ 
dures  are  robust/  resistant  methods  for  establishing  confi¬ 
dence  limits  on  means;  they  tend  to  down-weight  extreme 
observations  that  appears  as  outliers. 

For  small  sample  size,  (i.e.  n=15)  Winsorizing  at  g=2 
and  g=3  levels  provide  an  adequate  coverage.  3ut,  the  point 
availability  estimate  and  average  confidence  length  remain 
inadequate.  In  general,  the  Biweight  statistical  procedure 
caused  reduction  of  the  coverage  factor,  and  it  is  not  yet 
considered  effective  for  compensation. 

Next,  a  grouped  jackknife  procedure  was  applied  to  these 
cases. Simply,  up  and  down  times  data  were  grouped  two  by 
two,  and  average  of  these  groups  have  been  taken.  Obviously, 
this  process  reduced  the  sample  sizes  to  one  half  and  caused 
loss  of  degrees  of  freedom. However ,  smaller  variation  and 
more  stable  results  were  anticipated. Simulation  results  are 
presented  in  Tables  9  and  10.  Very  good  coverage  factors 
and  point  availability  estimates  ware  obtained,  but,  average 
confidence  length  remained  at  a  rather  high  level.  See 


r 


Appendix  C  for  histograms  of  both  cases  for  up  and  down 
times,  and  log- transform  of  down  times  and  pseudovalues 
after  grouped  jackknife  procedure  applied. 

Other  critical  simulation  results  were  obtained  from 
exponential  up  and  long-tailed  log-normal  h  down  times 
case.Por  both  JK  and  LNLJ  procedures,  above  statistical 
procedures  were  applied,  and  results  were  shewn  in  Tables 
11,  12,  13,  14,  15,  16. In  either  case,  after  implementation 
of  these  statistical  representations,  the  coverage  factor 
decreased.  Also,  see  Appendix  C  for  histograms  of  this 
case . 
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Table  1:  Simulation  Results  for  Distributional 

Sample  size=15 

Average  Variance 


Underlying 

Distribution 


A. 

Exponential 

up 

JK 

Log-normal  down 

LN 

MP 

B. 

Exponential 

up 

JK 

Exponential 

downLN 

MP 

C. 

Exponential 

up 

JK 

Gamma  down 

LN 

k=2 

MP 

D. 

Ditto 

JK 

k=  1/2 

LN 

MP 

E. 

Exponential 

up 

JK 

Long-tailed 

h 

LN 

h*0.  2 

MP 

?. 

Ditto 

JK 

h  =  0.  4 

LN 

MP 

Coverage 

Width 

Width 

0.  9380 

0. 0838 

0. 0022 

0. 9443 

0.0837 

0. 0016 

0.  96  98 

0.0959 

0.0026 

0.  9470 

0.0850 

0.0013 

0.  94  17 

0.1378 

0.0129 

0.  9868 

0. 1765 

0.0196 

0.  94  23 

0.0706 

0.0007 

0.94  19 

0.0764 

0.0010 

0.  9830 

0. 1001 

0.0031 

0.  9450 

0.1131 

0.0039 

0.  9161 

0.4574 

0.0932 

0.  9873 

0.5013 

0.0970 

0.  9256 

0.1171 

0.0097 

0.  9483 

0.1750 

0.0197 

0.  9780 

0.1951 

0. 0244 

0.  3753 

0.1663 

0.0343 

0.  9406 

0.2144 

0.0342 

0.9513 

0.2092 

0. 0321 

Assumptions 

Point 

Availability 

0. 9458 
0.9459 
0.9556 

0. 9471 
0.9312 
0.9421 

0.9465 
0. 9437 
0.9538 

0.9442 
0.8064 
0. 83  11 

0.9437 

0.9287 

0.9412 

0.9387 
0.9314 
0.  94  31 
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G. 

Long-tailed 

h 

upJK 

0.  9245 

0.  1077 

0.0024 

0.9428 

Exponential 

downLN 

0.  9058 

0.  1650 

0.0150 

0.9256 

h  =  0.  2 

UP 

0.  9803 

0.1973 

0. 0212 

0.  93  81 

H. 

Ditto 

JK 

0.  8837 

0.1470 

0.0048 

0.9351 

h=0.4 

LN 

0.  8552 

0.2100 

0.0179 

0.9163 

HP 

0.  9437 

0.2377 

0.0244 

0. 9284 

I. 

Exponential 

up 

JK 

0.  9476 

0.064  1 

0.0005 

0. 9483 

Weibull  down 

LN 

0.  9449 

0.0672 

0.0007 

0.9463 

k=2 

HP 

0.9838 

0.0917 

0.0023 

0.9501 

J. 

Ditto 

JK 

0.  9142 

0.1764 

0.0239 

0. 9404 

k=  1/2 

LN 

0.  9336 

0.5752 

0. 0921 

0.7649 

MP 

0.9751 

0.5319 

0.0966 

0.8464 

K. 

Exponential 

up 

JK 

0.8968 

0.1385 

0. 0220 

0. 9416 

Long -tailed 

lo 

g  ln 

0.  91  19 

0.1217 

0.0091 

0. 9483 

normal  h  down 

HP 

0.  91  83 

0.1113 

0.0058 

0.9576 

h-0.  2 

L. 

Ditto 

JK 

0.  9794 

0.1578 

0.0305 

0.9394 

h=0.  4 

LN 

0.  9009 

0.  1394 

0.01 36 

0.9493 

HP 

0.  9001 

0.1210 

0.0079 

0 . 9575 
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Table  2:  Simulation  Sesults  for  Distributional  Assumptions 


Sample  size=25 


Underlying 


Distributi 

on 

Cover  age 

Exponential 

up 

JK 

0.  9436 

Log- normal  down 

LN 

0.9473 

HP 

0.  9740 

Exponential 

up 

JK 

0.  9473 

Exponential 

downLN 

0.  9228 

HP 

0.9945 

Exponential 

up 

JK 

0.  94  75 

Gamma  down 

LS 

0.  94  14 

k=2 

HP 

0.  9895 

Ditto 

JK 

0.  94  64 

k  =  1/2 

LS 

0.8187 

HP 

0.  9884 

Exponential 

up 

JK 

0.  9313 

Long-tailed 

h 

LN 

0.  9464 

h=0.  2 

HP 

0.  9867 

Ditto 

JK 

0.  3784 

h  =  0.  4 

LN 

0.  94  75 

HP 

0.  9593 

Average 

Width 

Variance 

Width 

Point 

Availabili 

0. 0594 

0.0007 

0.9470 

0.0589 

0.0004 

0.9471 

0.0713 

0. 0008 

0. 9571 

0. 0597 

0.0004 

0.9471 

0.0978 

0.0052 

0.9322 

0.1302 

0.0086 

0.9450 

0.0506 

0.0002 

0. 9474 

0.0547 

0.0003 

0.9447 

0.0751 

0.0010 

0.9552 

0.0766 

0.0009 

0.9463 

0.4008 

0. 0790 

0.8168 

0.4606 

0. 0841 

0. 8435 

0.0834 

0. 0041 

0. 9456 

0.1213 

0.0  077 

0.9323 

0.1426 

0. 0108 

0.9450 

0.  1261 

0.0212 

0.9421 

0.1464 

0.0132 

0.9366 

0.1510 

0.0139 

0. 9481 
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Long- tailed 

h 

upJK 

0.9297 

0.0759 

0.0007 

0. 9454 

Sxponen  tial 

downLN 

0.  8901 

0.  1167 

0.0066 

0.9300 

h=0.  2 

HP 

0.98  34 

0.1447 

0.0099 

0.9422 

Ditto 

JK 

0.8802 

0. 1053 

0.0017 

0.9394 

h=0.  4 

LN 

0.8248 

0.  1521 

0.0079 

0.9228 

HP 

0. 9575 

0.  1771 

0. 0116 

0.9336 

Exponential 

up 

JK 

0.9448 

0.0464 

0.0002 

0. 94  89 

tfeibull  down 

LH 

0.94  3  3 

0.0485 

0.0002 

0. 9470 

k=2 

HP 

0.9879 

0.0691 

0.0008 

0.9494 

Ditto 

JK 

0.9222 

0.1188 

0. 0087 

0.9446 

k=  1/2 

LH 

0.8456 

0.5116 

0.0796 

0. 7765 

HP 

0.  97  27 

0.5461 

0.0854 

0.8245 

Exponential 

up 

JK 

0.  9010 

0. 1032 

0. 0123 

0.9440 

Long-tailed 

lo 

g  LN 

0.  91 41 

0. 0817 

0.0026 

0.9507 

noraal  h  down 

HP 

0.9078 

0.0792 

0.0018 

0.9600 

h  =  0.  2 

Ditto 

JK 

0.8835 

0.1184 

0.0183 

0.9426 

h  =0 . 4 

LN 

0.  9003 

0.0918 

0.0037 

0.9511 

HP 

0.8930 

0.0848 

0. 0024 

0.9604 
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Table  3 


:  Imp lamentation  of  "Winsorizing"  on  Exponential  Op 
and  Gamma  (k*1/2)  Down  Times  for  LNLJ  Procedure 
Sample  size*15 


Aver  age 

variance 

Point 

Coverage 

Width 

Width 

Availability 

G=  0 

0.9103 

0. 4212 

0. 0809 

0.8239 

G=  1 

0.9095 

0.3624 

0.0572 

0.8686 

G=  2 

0.9272 

0.3467 

0.0567 

0.8840 

G=  3 

0.9400 

0.3427 

0.0609 

0.8954 

Table  4:  Implementation  of  "Winsorizing"  on  Exponential  Op 
and  Gamma  (k»1/2)  Down  times  for  LNLJ  Procedure 
Sample  size=25 


Coverage 

Average 

Width 

Variance 

Width 

Point 

Availability 

G=  0 

0.8256 

0.3691 

0. 0670 

0.8319 

G=2 

0.8010 

0.2565 

0.0265 

0.8829 

G*3 

0.8361 

0.2456 

0.0247 

0.8908 

G*4 

0.8717 

0.2391 

0.0247 

0.8972 
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Table  5:  Implementation  of  "Winsorizing"  on  Exponential  Up 


and  Weibull 

( k=1/2) 

Down  Times  for 

LNLJ  Procedure 

Sample 

size=15 

Average 

Variance 

Point 

Coverage 

Width 

Width 

Availability 

G=0 

0. 9336 

0.5752 

0. 0921 

0.7649 

G  =  1 

0. 9418 

0.5056 

0.0828 

0. 8323 

G=2 

0.  9526 

0.4699 

0.0855 

0. 861 1 

G=3 

0.9510 

0.4464 

0.0907 

0.8836 

Table  6: 

Implementation  of  w 

Winsorizing"  on 

Exponential  Up 

and  Weibull 

(k=1/2) 

Down  Times  for 

LNLJ  Procedure 

Sam  pi 

e  size=25 

Average 

Variance 

Point 

Coverage 

Width 

Width 

Availability 

G=0 

0.  8456 

0.5116 

0.0796 

0.7765 

G=2 

0.  8692 

0.3874 

0.0516 

0.8537 

G=  3 

0.  9051 

0.3632 

0.0513 

0.8689 

G=4 

0. 9316 

0.3431 

0.0511 

0.8829 
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Table  7:  lap lementat ion  of  ••Biweighting"  on  Exponential  Op 
and  Samoa  (k=1/2)  Down  Times  for  LNLJ  Procedure 


Sample 

size 

Coverage 

Average 

width 

Variance 

width 

Point 

Availability 

Original 

1  5 

0.9103 

0.4212 

0.0309 

0. 8223 

Biweighted 

15 

0.3935 

0.2621 

0. 0448 

0.9117 

Original 

25 

0.8256 

0.3691 

0.0670 

0.8319 

Biweighted 

25 

0.3786 

0. 1768 

0.0159 

0.9195 

Table  8: 

Implementation 

of  ’'Biweighting"  on 

Exponential  Op 

and  Weibull  (k= 

1/2)  Down 

Times  for 

LNLJ  Procedure 

Sample 

Average 

Variance 

Point 

size 

Coverage 

width 

width 

Availability 

Original 

15 

0.9336 

0.5752 

0.0921 

0.7649 

Biweighted 

15 

0.8423 

0.3163 

0.0784 

0.9071 

Original 

25 

0.3456 

0.5116 

0.0796 

0.7765 

Biweighted 

25 

0.8374 

0.2134 

0.0362 

0. 9231 
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Table  9 


Table  9: 

Implementation  of 

"Grouped 

jackknife" 

Procedure  for 

Exponen 

tial  Op  and 

Weibull 

(k=1/2)  Down  Times 

Samp  le 

Average 

Variance 

Point 

Method 

Size 

Covera  ge 

Width 

Width 

Availability 

JK 

14 

0.9209 

0.  21  36 

0.  0353 

0.9381 

LS 

14 

0.9486 

0.3648 

0.0774 

0.9038 

21 P 

14 

0.9390 

0.3083 

0. 0672 

0.9073 

JK 

24 

0.9223 

0. 1319 

0.0118 

0.9429 

LS 

24 

0.9500 

0.2469 

0.  0400 

0.9148 

HP 

24 

0.9685 

0.25  18 

0.0452 

0.  9183 

Table  10:  I aplementati on  of  "Grouped  Jackknife"  Procedure 
for  Exponential  Op  and  Gamma  (k=1/2)  Down  Times 


Merhod 

Sample 

Siz  e 

Coverage 

Average 

Width 

Variance 

Width 

Point 

Availability 

JK 

14 

0.9426 

0.  1330 

0. 0080 

0. 9441 

LS 

14 

Q. 953  4 

0.2266 

0.0381 

0.9243 

HP 

14 

0.9608 

0.2186 

0.  0419 

0.  9262 

JK 

24 

0.9472 

0.0839 

0.0015 

0.9461 

LS 

24 

0.9462 

0.1448 

0.0174 

0.9293 

MP 

24 

0.9^83 

0.  1710 

0. 0259 

0.9322 
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Table  11:  Implementation  of  "Winsorizing"  on  Exponential  Op 
and  Long-Tailed  Log-Normal  (h*0.2)  Down  Times 
Sample  size=l5 


Coverage 

Average 

Width 

Variance 

Width 

Point 

Availability 

JK 

o 

ii 

u> 

0 .8968 

0.1385 

0.0220 

0.9416 

G=  1 

0.8190 

0.0700 

0.0019 

0.9544 

G=  2 

0.7450 

0.  0585 

0.001 1 

0.9560 

G=  3 

0.7160 

0.054 

0.0009 

0.9565 

LN 

G=0 

0.9119 

0.1217 

0.0091 

0.9483 

G=  1 

0.8560 

0. 0898 

0.0040 

0.9549 

G=2 

0.8040 

0.0759 

0.0029 

0.9575 

G=  3 

0.7460 

0.0668 

0.0027 

0.9598 
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12: 

I mplement 

ation  of  " w i 

nsorizing"  on 

Exponential  Op 

and  Long- 

Tailed  Log-Normal  (h=0.2) 

Down  Times 

Sample 

size=25 

Average 

Variance 

Point 

Coverage 

Width 

Width 

Availability 

JK 

G  =  0 

0.9010 

0.  1032 

0.0123 

0.9440 

G  =  2 

0.7330 

0. 0460 

0.0004 

0.9554 

G=3 

0.6950 

0.0416 

0.0001 

0.9563 

G=  4 

0.6620 

0. 0387 

0.0003 

0.9563 

LN 

G  =  0 

0.9141 

0. 0817 

0.0026 

0.  9507 

<3  =  2 

0.8240 

0.0594 

0.0010 

0.9553 

<3=3 

0.7750 

0.0539 

0.0008 

0.  9573 

G=4 

0.  7380 

0. 0488 

0.0006 

0.9589 
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Table  13:  Implementation  of  “Winsorizing"  on  Exponential  Up 
and  Long-Tailed  Log-Normal  (h=0.4)  Down  Times 
Saapla  size=l5 


Average 

Variance 

Point 

Coverage 

Width 

Width 

Availabil 

OK 

G=0 

0.8794 

0.1 578 

0.0305 

0.9394 

G=1 

0.7560 

0.0730 

0.0026 

0.9542 

G=  2 

0 .6870 

0.0588 

0.0015 

0.9560 

G=  3 

0.6460 

0.0537 

0.0012 

0.9565 

LN 

G=0 

0.9009 

0.1394 

0.0136 

0.9483 

G-  1 

0.8470 

0.  1049 

0.0069 

0.9543 

G=2 

0.7800 

0.0866 

0.0053 

0.9576 

G=  3 

0 .7180 

0.0763 

0.0048 

0.9603 
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Table  14:  I nplenentat ion  of  "Winsorizing"  on 
and  Long-Tailed  Log-Normal  (h~0.4) 


Coverage 

Sa  mple 

Average 

Width 

size=25 

Variance 

Width 

JK 

G=0 

0.8835 

0.1  184 

0.0183 

G=2 

0.6810 

0.0473 

0.0007 

G=  3 

0.6220 

0.0411 

0.0004 

G=4 

0.5750 

0.0380 

0.0004 

LU 

G=0 

0.9003 

0.0918 

0.0037 

G=2 

0.8000 

0.0666 

0.0018 

G=3 

0.7660 

0.0585 

0.0014 

G=4 

0.7000 

0.0  531 

0.0012 

Exponential  Op 
Down  Tines 

Point 

Availability 

0.9426 

0.9560 

0.9572 

0.9578 

0.9511 
0.9568 
0.9590 
C. 9609 
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Table  15:  Implementation  of  "Biweighting”  on  Exponential  Op 
and  Long-Tailed  Log-Normal  (h=0.2)  Down  Times 


Sa  mple 

size 

Cover  age 

Average 

width 

Variance 

width 

Point 

Availability 

JK 

Original 

15 

0. 8968 

0. 1385 

0.0220 

0.9416 

Biweighted 

15 

0.5870 

0.0528 

0.0010 

0.9554 

original 

25 

0.9010 

0. 1032 

0.0123 

0.9440 

Biweighted 

25 

0.6360 

0.0362 

0.0002 

0.9569 

LN 

Original 

15 

0. 91  19 

0. 1217 

0.0091 

0.9483 

Biweighted 

15 

0.70  10 

0.0604 

0.0019 

0.9607 

Original 

25 

0.9141 

0.0817 

0.0026 

0.9507 

Biweighted 

25 

0.6170 

0.0406 

0.0004 

0.9627 
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Table  16:  Implementation  of  '•Biweighting*’  on  Exponential  Op 


and  Long-Tailed 
Sample 

size  Coverage 

Log-Normal 

Average 

width 

(h=0.4) 

Variance 

width 

Down  Times 

Point 

Availability 

JK 

Original 

15 

0.  87  94 

0.1578 

0.0305 

0.9394 

Biweighted 

15 

0.6210 

0.0499 

0.0010 

0.9578 

Original 

25 

0.8835 

0.1184 

0.0183 

0.9426 

Biweighted 

25 

0.5500 

0.0340 

0.0002 

0.9590 

LN 

Original 

15 

0.  90  09 

0.1394 

0.0136 

0.9483 

Biweighted 

15 

0.6350 

0.0614 

0.0025 

0.9642 

Original 

25 

0.9003 

0.0918 

0.0037 

0.951  1 

Biweighted 

25 

0.5261 

0.0394 

0.0005 

0.9663 

IV.  CONCLUSIONS 


The  LMLJ  procedure  works  well  when  the  actual  down  times 
are  Exponential  or  Log-noraal.  However,  the  procedure  is 
very  sensitive  to  the  sample  variance.  If  the  sample  vari¬ 
ance  tends  to  vary  greatly,  sometimes  becoming  excessively 
large,  then  the  method  fends  to  fail;  the  intervals  shifted 
down.  For  instance,  long  tailed  Samma,  Weibull  and  Log-nor¬ 
mal  distributed  down  tines  provide  extremely  large  sample 
variances,  and  under  these  circumtances  the  method  fails. 

For  further  study,  one  might  apply  jackknifing  leaving 
out  one  observation  at  a  time,  instead  leaving  out  a  pair. 
Obviously,  this  is  going  to  increase  the  degrees  of  freedom, 
but  might  by  help  reducing  the  variance  and  the  confidence 
interval  length.  On  the  other  hand,  grouping  more  than  two, 
might  provide  more  stable  results. 

In  addition,  in  order  to  see  what  happens  in  one-sided 
availability  estimation.  Exponential  up  and  Log-normal  down 
and  Exponential  up  aeibull  (k=1/2)  down  times  were  simulated 
and  results  are  shown  in  Appendix  D. 


APPENDIX  A 


COMPOTES  PROS  RAMS 

Simulation  program  consists  of  main  program  and  several 
subroutines.  Main  program,  for  JK,  LNLJ  and  MP  procedures 
computes  availability  confidence  limits  and  point  availabil¬ 
ity  and,  scores  the  coverage  for  each  replication.  Then, 
after  1000  replications  computes  the  statistics  of  these 
parameters  and  prints  out  the  results  for  given  underlying 
distributions. 

Subroutine  CONP  computes  two-sided  confidence  limits  of 
a  given  data  vector.  Subroutine  MOMENT  computes  the  fourth 
moment  of  given  data  vector,  subroutine  LNLJ  generates  the 
up  and  down  times  from  given  underlying  distributions,  and 
takes  the  log  transform  of  these,  and  computes  the  means  ana 
the  variances  of  these  data  vectors. 

Subroutine  BIWGT  computes  the  Biweight  estimates  of  a 
given  data  vector,  and  uses  Subroutine  MEDIAN  fcr  computing 
the  median  of  this  data.  Subroutine  WINSOR  wir.sorize  a 
given  data  vector  for  given  "g"  level. 
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APPENDIX  3 


LONGlSAILSD  LCS- normal  distribution 

(1)  Let  H  be  log-normal  random  variable  which 
lnW^N(v#a2).  k  moment  of  M  as  follows, 

EC  W k  3=  exp  (k  u  *  1/2*  k2  j2), 

see  [Ref.  9]. 

So, 

E[W  ]=  exp  (  u  ♦  1/2*  a2) 

E[ '4 2  ]■  exp  (2  u  +  2  a2) 

(2)  Create  a  stretched  log-normal  random  variable 
D»  c«ff*(1*h«V) 

where  c  and  h  are  constant  and  0<h<l. 

E[D  ]=  c(E[  W  ]  +  h*E[  w2  ]} 

E[D]=  c  (exp  (u+ 1/2»  o  2)  ♦  h«exp(2,j  *2  a2)} 

c=  E[D]  (  exp  (  u  +1/2*  a2)  ♦h*exp(2u+2a2)}-1  (a .  1 ) 

(3)  On  the  other  and,  we  can  write  the  median  for  which 
is  the  function  of  the  expected  value  as  follows, 

H*  d  (0.5)  =  V  S[  D] 

a*  c«exp  (  u)  •  { 1+h*exp  (  y))=  *E[D]  (a.  2) 
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where  ♦  is  the  median  location  factor  of  the  new  dis¬ 
tribution.  If  we  put  equation  (a .  1 )  in  equation 
(a.  2)  ,  we  obtain 


>l>  •  21  D  ]=  E[  D  ] 


exp  (u)  •  (1  +  h«exp(|i)) 


exp  (  u+1/2  •  a2)  ♦  h«exp(2y*2  c2) 

After  the  cancellations,  we  obtain 
!f/(exp(u+  1/2*  a  2)  ♦  h»exp(2y  ♦  2  a2)  }=  1  +  h*exp  (y ) 
^  «exp  (1/2  3  2)  -  1  =  h*exp  (y)  •  (1  -  i|i  *exp  (2  a  2)  ) 

1- *  *exp (1/2 •  3  2) 
h=  exp(-y) - - 


■4*  *exp  (2  j  2)  -  1 
If  h=0,  then  ii  -  exp  (-1/2*  32) 

1  ♦  exp  (y) 

if  h=  1,  then  ip  3  -  ■  - - - 


9XP(U*2  J2)  *■  exp  (1/2»  a2) 

So,  if  we  get 

exp(-1/2*a2)  <  <  exp(y*2  a2)  *exp  (1/2*  a2) 

then  we  obtain  h  between  0  and  1. 
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A££ENDI2  C 


DISTRIBUTIONAL  ?I22BI5 

The  LHLJ  method  fails  for  Exponential  up  and  Gamma 
(k=1/2)  down.  Exponential  up  and  Weibull  (k=1/2)  down,  and 
Exponential  up  and  Long-tailed  Log-  normal  h  down  times 
cases.  In  order  to  demonsrate  the  distributional  charecter- 
istics  of  this  cases,  histograms  are  obtained  from  simula¬ 
tion  for  up  times,  down  times,  log  transformations  of  down 
times  and  pseudovalues.  Also,  histograms  are  repeated  after 
grouping  the  observations  in  order  to  observe  the  grouping 
effect  on  the  parameters.  Distributional  charecteristics  of 
these  parameters  such  as  means,  variances,  moments  and 
shapes  can  be  observed  from  these  histograms. 
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(Jamma  (k  =  1/2)  down  times  after  grouping. 
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Long-tailed  Log-normal  h  down  times 
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mg-tailed  log-normal  down  times  after 


FPECUffcCIES  SAMPLE  SUE  -  $00 


rw*  ■  ■■ 


FPEOOiMCIES  SAMPLE  SUE  >  500 


21.  Pseudovalues  of  JK  procedure  for  Exponential  up  and  long¬ 
tailed  log-normal  down  times  after  grouping. 


APPENDIX  D 


SllfcSIfiED  hUlU&lltlll  SHIHAIIQS 
Best  and  worst  cases  for  LNLJ  procedure  were  simulated 
for  one-sided  availability  estimation,  and  results  are  shown 
in  Tables  17  and  18  as  follows; 


Table 

17:  Exponential  Op 

and  Log- 

normal 

Down  Ti 

mes 

N 

=  15 

N 

=  25 

JK 

LN 

JK 

LN 

Low. 

Availability  Bound 

Cover 

0.999 

1.000 

0.999 

1.000 

Opp. 

Availability  Bound 

Cover 

0.933 

0.899 

0.966 

0.965 

Avg. 

Lower 

Availability 

3ound 

0.889 

0.879 

0.910 

0.  908 

Avg. 

Upper 

Availability 

Sound 

0.966 

0.964 

0.966 

0.966 

Table  13:  Exponential  'Jp  and  Weibull  (k=1/2)  Down  Times 

N= 1 5  N=  25 


JK 

LN 

JK 

LN 

Low.  Availability  Bound 

Cover 

0.  300 

0.9  90 

0.9  14 

0.997 

Opp.  Avail 

ability  Bound 

Cover 

0.966 

0.853 

0.955 

0.705 

Avg.  Lower 

Availability 

Bound 

0.344 

0.460 

0.331 

0.  500 

Avg.  Upper 

Availability 

Bound 

0.972 

0.970 

0.975 

0.955 
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Notica  the  low  Average  Lower  Availability  Sound  provided 


by  the  LNLJ  procedure  for  the  Weibull  (lc*1/2),  Sesult 
that  the  one-sided  availability  estimation  behaves  as 
two-sided  availability  estimation. 
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